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ABSTRACT 

Context. Classical T Tauri stars are young solar-type stars accreting material from their circumstellar disks. Thanks to a favorable 
inclination of the system, the classical T Tauri star AA Tau exhibits periodic optical eclipses as the warped inner disk edge occults the 
stellar photosphere. 

Aims. We intend to observe the X-ray and UV emission of AATau during the optical eclipses with the aim to localize these emitting 
regions on the star. 

Methods. AATau was observed for about 5h per XMM-Newton orbit (2 days) over 8 successive orbits, which covers two optical 
eclipse periods (8.22 days). The XMM-Newton optical/UV monitor simultaneously provided UV photometry (UVW2 filter at 206 nm) 
with a ~15min sampling rate. Some V-band photometry was also obtained from the ground during this period in order to determine 
the dates of the eclipses. 

Results. Two X-ray and UV measurements were secured close to the center of the eclipse (AV ~ 1.5 mag). The UV flux is the highest 
just before the eclipse starts and the lowest towards the end of it. UV flux variations amount to a few 0. 1 mag on a few hours timescale, 
and up to 1 mag on a week timescale, none of which are correlated with the X-ray flux. We model it with a weekly modulation (inner 
disk eclipse), plus a daily modulation, which suggests a non-steady accretion, but needs a longer observation to be confirmed. No such 
eclipses are detected in X-rays. Within each 5 h-long observations, AA Tau has a nearly constant X-ray count rate. On a timescale 
of days to weeks, the X-ray flux varies by a factor of 2-8, except for one measurement where the X-ray count rate was nearly 50 
times stronger than the minimum observed level even though photoelectric absorption was the highest at this phase, and the plasma 
temperature reached 60 MK, i.e. a factor of 2-3 higher than in the other observations. This X-ray event, observed close to the center 
of the optical eclipse, is interpreted as an X-ray flare. 

Conclusions. We identify the variable column density with the low-density accretion funnel flows blanketing the magnetosphere. 
The lack of X-ray eclipses indicates that X-ray emitting regions are located at high latitudes. Furthermore, the occurrence of a strong 
X-ray flare near the center of the optical eclipse suggests that the magnetically active areas are closely associated with the base of the 
high-density accretion funnel flow. We speculate that the impact of this free falling accretion flow onto the strong magnetic field of 
the stellar corona may boost the X-ray emission. 

Key words. Stars: individual: AATau - Stars: pre-main sequence - Stars: flare ~ X-rays: stars - Accretion, accretion disks 

1 . Introduction the X-ray luminosity is proportional to the plasma emission mea- 
sure, which scales linearly with the plasma volume but with the 

T Tauri stars (TTSs), i.e. young (1-lOMyrs) solar-type stars, gqu^re of the plasma electronic density. That most of TTSs X- 

are conspicuous X-ray emitters. Their high X-ray lummosities j-ay emission arises in an active magnetic corona is supported 

(Lx - 10 ergs ) compared to the Sun (Lx - 10 ergs by direct Zeeman measurements on photospheric spectral lines 

at solar maximum), and intense flaring activity (up to Lx ^ which indicate surface magnetic fields of a few kilogauss (e.g., 

IQJz-JJ erg s-i) make them appear as extremefly active young Quenther elal. 19 99; Johns-Krull eLal. 1999, 2004; Johns-Krull 

Suns m the X-ray domain. The analogy with the solar activ- 2007; Y an7^a[Mw~± However, the dynamo mechanism pro- 

ity has been quite successful in ascribing the X-ray emission of Cueing the ma gnetic field in fli ese f ully convectiv e stars is still 

TTSs to an optically thin, magnetically confined coronal plasma discussed (e.g.. lPreibisch itd]l2CioaiBriggs et al.ll2007l) . 
in collisional equilibrium at temperatures of 10-100 MK, emit- 
ting a thermal bremss trahlung continuum and emission lines The solar paradigm, where the X-ray emitting plasma is con- 

(see, e.g., review by iFeigelson & Montmerl ^ 119991). As ob- fined in magnetic loop with both feet anchored on the steflar pho 



serve d in the X-ray coronae of active stars (e.g.. iNess et all tosphere, has been questionned in the context of a Classical TTS 



|2004|), the enhanced X-ray luminosity of TTSs can easily be ex- (CTTS), which accretes material from its circumstellar disk, 
plained by coronal structures with high plasma density, because Inside a few stellar radii above the CTTS surface, the stellar 
magnetic field pressure is larger than the ram pressure of the ac- 



* Figures|2l[3]and[9l and Appendix |A| are only available in electronic creting gas. As a result, the steUar magnetosphere truncates the 
form via http : //www . edpsciences . org . inner accretion disk and controls the accretion flows. The gas is 
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mainly accreted from the disk edge to the stellar surface along 
the dominant large scale stellar magnetic lines, creating accre- 
tion funnel flows. The free-falling gas hits the stellar surface at 
the feet of the accretion funnel flows, where the kinetic energy is 
dissipated in a shock producing hot excess emission (s ee review 
on magnetospheric accretion by iBouvier et al ]|2007bl) . 

The X-ray grating spectrometers aboard Chandra and XMM- 
Newton are able to obtain spectra of the X-ray brightest CTTSs, 
where the emission line triplets of He-like elements are re- 
solved, which provides a powerful tool to ass ess the elec- 
tronic density of the X-ray emitting plasma ("e.g.. iPorquet et all 
l200lb . In several CTTSs, plasma with high electronic density 



Table 1. Journal of the XMM-Newton observations of 
AATau (PI: J. Bouvier). 



("e 



10 cm ) and low temperature (~3MK), untypical 



of stellar coronae , were identified, and therefore attributed to 
accretion shocks (Kastner et al."2002'; 'Stelzer & Schmitt"20()4t 
Schmittetal. 2005; Robrade & Schmitt 2006; Argiroffiet aL| 
20071: IGiinther et alj 120071: ITelleschi et all 120071)"; whereas 



some CTTSs display no such evid ence jAudard et alj 120051; 
ISmith et al.ll2005l;lGiidel etaHIIOOTbl) . 

During the Chandra Orion Ultradeep Project (COUP, see 

iGet man et al. 2005b), where the TTSs of the Orion nebula clus- 
ter were monitored nearly continuously over ~13 days with the 
Advanced CCD Imaging Spectrometer, numerous X-ray flares 
were observed. In a few cases, a size of several stellar radii 
was derived for the magnetic loop confining the X-ray emitting 
plasma, large enough to connect the stel lar surface with the edge 
of the accretion disk dFavata et al.ll2005h . 

We propose to directly constrain the source and location of 
the X-ray emission in CTTSs by using eclipses. Eclipse map- 
ping have been successf ully used in binary activ e stars to recon- 
struct the coronae (e.g.. lGtidel et alj 200 iL 2003,) , or to localize 



the fl aring plasma (e.g., I Schmitt & Favata 



T99I ISchmittltaP 



I2003h . Our target star is AA Tau, located in the Taurus molec- 
ular c loud complex at a distance of ~140pc (e.g.. lKenvon et all 
11994 *). AA Tau is a quite typical member of the CTTS class, with 
a K7 spectral type, a bolometric luminosity of ~O.8L0, a stel- 
lar mass of ~O.8M0, and a stellar radius of ~1.85Ro; exhibit- 
ing moderate accretion disk diagnostics (near-IR excess, optical 
veiling, Balmer line emission), but_with the rema rkable prop- 
erty to be viewed nearly edge-on. iBouvier et alj (1 19991 12003L 
l2007ah reported evidence for a modulation of the photospheric 
flux and spectroscopic diagnostics with a period of 8.22 days, 
corresponding to the rotational period of the star This stellar 
flux modulation was interpreted as the periodic eclipse of the 
stellar photosphere by the optically thick, magnetically-warp in- 
ner disk edge located at 8.8 stellar radii. Photopolarimetric vari- 
ations confirm the presence of an optically thick wall located at 
the disk edge, and e chpsing periodically the stellar photosphere 
dMenard et alj|2003h . 

We obtained 8 obser vations of AATau with XMM- 
Newton ("Jan sen et al.ll200lh. which a llows simultaneous obser- 
vations with an EPIC pn dStriider et al. 2001) and two EPIC 
MOS ( Turne ret al.l 1200 iT) X-ray spectroimaging ca meras, and 
the Optical/UV monitor (OM; iMason et all 1200 ll) . We sup- 
plemented these X-rays and UV observations with an optical 
ground-based monitoring of AA Tau to secure the dates of the 
optical eclipses. Th ese XMM- Newton obser vations were previ- 
ously reported in Schmitt & Robradd (l2007l hereafter SR), who 
used the minimum of the UV light curve as proxy of the optical 
eclipse. SR found variable X-ray absorption "such that the times 
of maximal X-ray absorption and UV extinction coincide". SR 
introduced an additional absorption in a disk wind, or a peculiar 
dust grain distribution to reconcile the high value of the X-ray 
absorption outside the eclipse and the low optical extinction. 



Obs. 


Rev. 


Obsid 


Feb. 2003 


Exposure" 








W 


W 


1 


583 


015680201 


14.10-14.29 


4.7 


2 


584 


015680301 


16.13-16.33 


4.7 


3 


585 


015680501 


18.07-18.26 


4.7 


4 


586 


015680401 


20.02-20.22 


4.7 


5 


587 


015680601 


22.13-22.35 


54 


6 


588 


015680701 


24.03-24.25 


5.3 


7 


589 


015680801 


26.52-26.71 


4.6 


8 


590 


015680901 


28.30-28.50 


4.7 



The observation beginning, end, and duration is given for MOS 1. 



In Sect. |2] we present the X-ray and UV properties of 
AA Tau based on a reanalysing of the full data set provided by 
these XMM-Newton observations. In particular, we report X-ray 
(pn+MOS 1 +MOS2) and UV light curves with a time resolution 
of ~15min. We show that a bright and hot flare was observed 
during the second observation. In Sect. [3] thanks to our ground- 
based observations and optical OM data, we determine the dates 
of the optical eclipses, which allow us to compare the UV and X- 
ray variabilities versus the rotational phase. We show that the UV 
minimum is outside the eclipse, and that the (confirmed) varia- 
tion of the column density are not correlated with the rotational 
phase. In Sect.|4] we propose another origin for this variable col- 
umn density, and we discuss the origin of the X-ray flare that was 
observed during an optical eclipse. 

2. X-ray and UV properties of AATau 

2.1. XMM-Newton observations and event selections 

Our observational strategy with XMM-Newton was to cover two 
consecutive modulation periods (~17 days) in order to demon- 
strate the reproducibility of the phenomenon from one rotational 
cycle to the next. The temporal sampling of the X-ray light 
curve needs not to be very dense because AA Tau spends nearly 
the same amount of time being eclipsed as being entirely vis- 
ible. Therefore, we requested one 4h-exposure with EPIC pn 
per XMM-Newton orbit (2 days) over 8 successive orbits. We 
used the full frame science mode of the EPIC cameras with 
the medium optical blocking filter The pointing nominal co- 
ordinates were 04''34'"55!5, 24°28'54'.'0 (J2000 equinox). The 
journal of the XMM-Newton observations of AA Tau is given in 
TablelU 

The data reduction was made using the XMM- 
Newton Science Analysing System (SAS, version 7.0.0) 
For each observations, the event lists for each camera of EPIC 
were produced using the SAS tasks epchain and emchain, 
respectively. We used the background lightcurves computed by 
these task in the 7.0-15.0keV energy range to determine the 
time intervals affected by background proton flares. More than 
half of the observing time is affected by bad space weather. 
For each instrument, we made from the low background time 
intervals a sky image with 3"-pixels in the 0.5-7.3 keV energy 
rang4B The X-ray counterpart of AA Tau was detected in all the 
observations. 



' We selected single, double, triple, and quadruple pixel events (i.e. 
PATTERN in the to 12 range), and also applied the predefined filter 
#XMMEA_EM for MOS. 



#1 



#2 
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Day in February 2003 

Fig. 1. XMM-Newton observations of AA Tau. Top and bottom panels show the background subtracted UV ( 1 80-250 nm) and X-ray 
(0.5-7.3 keV) hght curves obtained with the Optical/UV Monitor (OM) and EPIC (pn+MOSl+MOS2), respectively. In each panel, 
the label indicates the time interval used to bin the light curve. Grey hourglasses and crosses data points indicate UV photometry 
obtained with the large and small OM central imaging window. The X-ray count rates are corrected from (circular) aperture. The 
gray stripes indicate the beginning and the end of the EPIC observation. 



For the X-ray light curves of AATau, we selected the 
sourceH-background events within a circular region centered on 
AA Tau . The extraction position and radius were optimized to 
maximize the signal-to-noise ratio in the sky image. The back- 
ground in MOS and pn was extracted using an annular region 
centered on AATau and a box region at the same distance to 
the CCD readout node, respectively, where areas illuminated by 
other weak X-ray sources were excluded. For the X-ray spec- 
tra, we used the same extraction regions and only time inter- 
vals with low background; we selected events with energy above 
0.3 keV and the usual stronger selection criteri^ We computed 
the corresponding redistribution matrix files and ancillary re- 
sponse files. 

2.2. X-ray light curves 

For each instrument and observations, we first built the 
sourceH-background and the background light curves with 1 s 
time bins starting at the first good time interval (GTI) of MOS 1 . 
We rebinned the light curve to 900 s to increase the signal. Then, 
we subtracted from the source+background light curve the back- 
ground light curve scaled to the same source extraction area. 

From the GTI extension, we scaled up with an IDL routine 
count rates and errors affected by any lost of observing time, 
mainly due to the triggering of counting mode during high flar- 
ing background periods, when the count rate exceeded the de- 
tector telemetry limit. We also corrected count rates and eiTors 
for circular aperture photometry using the fraction of PSF counts 
inside the (circular) extraction region calculated by the SAS as- 
suming a fixed photon energy of 1 .5 keV. Finally, the light curves 
of the three detectors were summed to produce the EPIC light 
curves. We estimated the missing pn data at the beginning of the 
observations by multiplying the MOS1H-MOS2 count rates by 



^ For pn, we selected only single and double pixel events (i.e. 
PATTERN in the to 4 range) with FLAG value equal to zero; for MOS, 
we selected PATTERN in the to 12 range, and applied the predefined 
filter #XMMEA.SM. 



1.17, the median scaling factor between MOS 1 -hMOS and pn in 
the second observation where AA Tau was the brightest. 

The bottom panel of Fig. [T] shows the EPIC light curves of 
AATau in the 0.5-7.3 keV energy range. Within each 5h-long 
observations, AATau exhibited a nearly constant X-ray flux. 
The minimum X-ray flux, that we will call hereafter the qui- 
escent level, was observed during the observation #5, where 
the averaged EPIC count rate was 0.030 + 0.002 counts s"'. 
On a timescale of days to weeks, the X-ray flux varies by a 
factor of 2-8, except between the end of the first observation 
and the beginning of the second observation, where the EPIC 
count rate jumped in less than 2 days from 0.093 + 0.003 to 
1.42 + 0.01 counts s"\ i.e. a level 47 + 3 times stronger than the 
quiescent level. Then, the EPIC count rate decayed in less than 2 
days to 0.075 + 0.003 counts s"', i.e. a level only 2.5 ± 0.2 times 
stronger than the quiescent level. 

Such large amplitudes in the X-ray fluxes of young stellar 
objects are usually observed during X-ray flares, which have 
typical light curves with fast rise and peak phase, and slower 
(exponential) decay phase, associated with fast h eating and slow 
coolin g of the magneticall y confined plasma (e.g.. llmanishi et al] 
20031; iFavata et all l2005h . X-ray flares w ith unusuaHy long 
rise phases have also been reporte d (e.g., iGrosso et alJ 120041 ; 
IWang et al.ll2007l;lBroos et alj|2007h . 

For a comparison purpose, we note that, assuming a 
typical convertion ratio betwee n XMM -Newton/EPlC and 
Chandra/ ACIS-I count rates (e.g.. lOzawa e t al. 2005), the aver- 
age EPIC count rate of AA Tau, ~0.08 counts s"' at the distance 
of the Taurus molecular cloud (d ~ 140 pc), would convert to 
~0.001 Qiandra/ACIS-I counts s ' at the distance of the Orion 
nebula (d ~ 450 pc). AATau, put in the Orion nebula cluster, 
would have then been brighter than 67% of the X-ray sources 
in the COUP. To check whether the behaviour of the light curve 
of AA Tau is consistent with the one of an X-ray flare, we can 
compare it to the light curve data set obtained in the COUP. 

We use the COUP results of the Bayesian block (BE) vari- 
ability analysis (developed by Scargle 1998; adapted for the 
COUP data set and coded in IDL by one of us, N.G.), which 
segmented the X-ray light curves into a contiguous sequences 
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Table 2. Best parameters of simultaneous fitting of EPIC pn, MOSl, MOS2 spectra with XSPEC. 



Plasma temperature Emission measure Goodness-of-fit Flux'' Intrinsic luminosity t/^ 



Obs. 


Rate" 


Nh 


Ti T2 


EMi EMi 




<3 




0.5-2, 


2-8, 0.5- 


~8keV 






(cuts s"') 


(10^2 cm-2) 


(MK) 


(10^3 cm-3) 




(%) 




(lO^^ergs- 


') 




1 


0.057 


0.95 


25.9 


1.0 


0.93 (50) 


61 


1.9 


0.6 


0.4 


0.9 


-3.5 




+0.002 


0.87-1.02 


23.9-28.4 


0.9-1.0 
















2 


0.800 


1.54 


43.9 


14.0 


1 f\r\ / c r\f\\ 

1.09 (599) 


5 


37.8 


8.2 


9.2 


17.4 


-2.2 




+0.008 


1.51-1.57 


42.7^5.2 


13.6-14.3 
















2 


0.800 


1.80 


14.6 59.4 


9.8 9.3 


1.03 (597) 


28 


37.9 


11.8 


9.5 


21.3 


-2.2 




±0.008 


1.75-1.85 


13.4-15.7 54.6-66.0 


8.8-10.7 7.8-10.8 
















3 


0.046 


1.36 


20.8 


1.1 


1.21 (30) 


20 


1.6 


0.7 


0.3 


1.0 


-3.5 




+0.002 


1.23-1.51 


18.4-23.4 


1.0-1.3 
















4 


0.031 


1.33 


21.5 


0.8 


0.79 (25) 


76 


1.1 


0.5 


0.2 


0.7 


-3.6 




±0.002 


1.17-1.52 


18.7-24.8 


0.6-0.9 
















5 


0.019 


1.04 


22.6 


0.4 


1.32(15) 


18 


0.7 


0.2 


0.1 


0.4 


-3.9 




±0.002 


0.85-1.27 


18.3-28.8 


0.3-0.5 
















6 


0.050 


1.54 


40.5 


0.9 


1.20 (39) 


18 


2.3 


0.5 


0.5 


1.1 


-3.5 




±0.002 


1.40-1.72 


34.6-47.4 


0.8-1.0 
















7 


0.147 


1.02 


44.8 


2.1 


1.05 (110) 


33 


6.5 


1.2 


1.4 


2.6 


-3.1 




±0.005 


0.97-1.08 


41.3^8.6 


2.0-2.2 
















8 


0.043 


1.53 


23.3 


1.2 


1.19(29) 


23 


1.8 


0.7 


0.4 


1.1 


-3.4 




±0.002 


1.39-1.69 


20.7-26.5 


1.0-1.4 

















Notes: We fit the X-ray spectra with the continuum and emission lines produced by an optically thin plasma in thermal collisional ionization 
equilibrium model (vapec). We use for the plasma element abundances typical values observed in the coronae of young stars with fine 
X-ray spectros copy (Gudeletal. 2007a, see Table \AA\. The wa bs photoelectric absorption model use photoionization cross sections of 
|m orrison & McCammon ( 198 3), and solar abundances of lAnders & Grevesse ( 198 9). Errors are given at the 68% confidence level (i.e. Ax^ = 1.0 
for each parameter of interest), that corresponds to Itr for Gaussian statistics. For observation #2, a better fit is obtained using a plasma with two 
temperature components. 

" pn count rate (0.2-12 keV). 
Observed X-ray flux (0.5-8.0 keV) in unit of 10~" erg cm'^ s"' . 

Logarithm of the X-ray intrinsic luminosity (0.5-8 keV) to the bolometric luminosity (0.8 Lq) ratio. 



of constant count rates (see iGetman et al.ll2005bl) . We define 
a source to be variable, if there is more than one BB; with 
BBjnin and BB^^^, the minim um and the rnaximu m count rate 
levels, respectively (see, e.g., IStassun et al.ll2007l for an appli- 
cation of COUP time-averaged X-ray variability). The latter and 
the former are viewed as the quiescent level and the peak level 
of the brightest flare, respectively. Applying this criteria, there 
are 977 variable sources (out of 1616 COUP sources). We find 
only 20 variable sources with peak amplitude (BBmax/Bfimin) 
and duration larger than 45 and 4.7 h., respectively, as ob- 
served for AA Tau. For comparison, our subsample sources have 
fifimax - 10"'*-2 counts s"', i.e., they are at peak between 200 
times fainter and 100 times brighter than AA Tau (at the distance 
of the Orion nebula) at peak. 

Then, we make a visual examination of the selected BB light 
curves to eliminate sources with peak flare BB having a spuri- 
ous long duration (including for example a passage through the 
van Allen belts), and/or sources with a decay phase too slow to 
reproduce the level observed at the beginning of the third ob- 
servation of AA Tau. For an exponential decay phase with a de- 
cay timescale tj, this latter criteria is equivalent to < 0.6 day. 
Finally, we exhibit the brightest flares from COUP 313, 874, and 
970, which fullfill our criteria (see online Fig.|2]l. These sources 
are bon a fine members of the Orion nebula cluster (iGetman et al.l 
l2005ah . 

We conclude that a bright X-ray flare with a rapid cooling 
phase can well reproduce the large amplitude, and also the flat- 
ness of the light curve of AA Tau at its maximum. The following 
X-ray spectra analysis confirms this interpretation. 



2.3. X-ray spectra and plasma parameters 

For each observations the pn, MOSl, and MOS2 spectra were 
binned to 25 counts per spectral bin. X-rays are detected up to 
lOkeV (see online Fig.O. The spectra are featureless, except in 
the second observation where a prominent line around 6.7 keV 
is visible both in the pn and MOS spectra, corresponding to the 
Fe XXV triplet emission line (Fig.|4|i. 

The pn, MOSl, and MOS2 s pectra were fitted simulta ne- 
ously with XSPEC (version 12.3.0: iDorman & Arnaudll200l[) to 
derive the plasma parameters. The model is an X-ray emission 
spectrum from collisionally-ionized diffuse gas, output from 
the Astrophysical Plasma Emission Code (vapecfl), that in- 
cludes continuum and emission lines. The plasma element abun- 
dances were fixed to typical values measured in the coronae 
of young stars with grating X-ray spectroscopy (see online ap- 
pendix Table lA.ll . to allow a direct comparison with the XMM- 
Newton E xtended Survey of the Taurus molecular cloud (XEST; 
iGiidel et al. 2007a) . This emission model is combined with the 
wabs photoelectric absorption model, which is based on the pho - 
toionization cross sections of 'Morrison & McCammorj ( Il983h . 
and the solar abundances of Anders & Grevesse ( 19891). The on- 
line Fig. [3] shows our best fits. For the second observation, we 
find a better fit by adding a second temperature component, 
which reduced the;\f^ from 655.4 (for 599 degrees of freedom; 
d.o.f.) to 617.2 (for 597 d.o.f.). Given the new and old values of 
and number of degrees of freedom, a F-test indicates a prob- 
ability of ~I0"^ for the null hypothesis; therefore, we conclude 

^ More information can be found at: 

|http://hea-www.harvard.edu/APEC/soiirces-apec.htnil| . 
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Fig, 4. EPIC pn (black), MOSl (red), MOS2 (green) spectra of 
AATau for observation #2. The lines show our best fit using 
two temperature plasma combined with photoelectric absorption 
(Table|2]i. The residuals are plotted in sigma units with error bars 
of size one. 



that it is reasonnable to add this extra temperature component 
to improve the fit. Our best fit with a two-temperature plasma is 
shown in Fig.H] 

Table |2] gives the corresponding plasma parameters. On 
timescale of 2-days, the photoelectric absorption of the X-ray 
spectra is not constant, but the observed relative variations are 
lower than a factor of two. The minimum and maximum values 
of the corresponding column density are ~1.0 x 10^^ cm"^ and 
1.8 X 10^^ cm"^, respectively; the latter was observed during the 
second observatiorQ. SR didn't specify the model of X-ray ab- 
sorption that they used, but their one-temperature plasma model 
gave qualitatively similar results. We stress that the value of the 
column density can even be increased by 50% when revised so- 
lar abundances (i.e., metal poor) are adopted for the abundances 
of the absorbing material (see online AppendixlAll. 

The plasma temperature is ~23MK during the low activ- 
ity levels. The observations showing an increase of X-ray count 
rates (namely #2, #6, and #7) also correspond to phases with 
the highest plasma temperatures (60 MK, 40 MK, and 45 MK, 
respectively), and therefore to flaring activity. 

During the second observation the X-ray surface flux, i.e., 
the intrinsic X-ray luminosity divided by the stellar surface, 
peaked to 10^ erg s"' cm"^. The temperatures of the hot (60 MK) 
and cool (15 MK) plasma components, associated with this ele- 
vated level of X-ray surface flux, are consistent with the ones 
observed in the mos t active TTSs of COUP (see Fig. 11 of 
iPreibisch et al.l l2005h . The cool plasma component is usually 
observed in the coronae of active stars, and may define a fun- 
damental coronal structure, which is probably r elated to a class 
of co mpact loops with high plasma density (see IPreibisch et al] 
I2005L and references therein). 

We make an estim ate of the hot loop length using the method 
of iReale et al.l d 19971) (see alsojFavata et al. 2005, for an appli- 
cation on COUP data). Assuming a peak temperature of 60 MK 
and a flare decay time lower than 0.6 day (see Sect. 12.21 ), we find 
for the semi-circular loop a half-length lower than 7.1 R* for a 
freely decaying loop, with no heating; or lower than 1.6 R* for 



The increase of the spectrum slope between 1 and 2 keV observed 
between observations #1 and #2, is not only due to the increase of the 
absorption, as argued in SR, but also for 20% to the increase of the 
plasma temperature. 



a strongly sustained heating in the cooling phase. Therefore, 
the height of the semi-circular loop is lower than 1^.5 R*. 
Therefore, this loop cannot connect the stellar surface and the 
inner accretion disk, distant by 7.8 stellar radii, and is likely an- 
chored on the stellar photosphere. 

The high temperature and X-ray surface flux observed dur- 
ing the second observation point to an enhanced X-ray activity 
produced by a bright X-ray flare. Such bright flares on the Sun 
are sometimes associated with coronal mass ejection (CME). 
Therefore, we cannot rule out that the maximum of column den- 
sity observed during the second observation is due to this ener- 
getic event. 

2.4. Comparison with previous X-ray observations of AA Tau 

AA Tau was previously observed several times in X-rays at dif- 
ferent epochs: on March 4, 1980 and February 7, 1981 with 
the IPC on board Einstein, with an exposure of 0.6 h and 2.8 h, 
respectively; on August 10, 1990 during the ROSAT All Sky 
Survey (RASS) with the PSPC on board ROSAT, with an expo- 
sure of 0.2 h; and on February 22 and August 18, 1993 during 
two pointed PSPC observations of the young binary Haro 6-13 
(PI: H. Zinnecker), with an exposure of 0.6 h and 1 .5 h, respec- 
tively. 

IWalter & Kuhil (Il98ll) reported a detection with 0.030 + 
0.005 IPC counts s ' , during the first observation with Einstein 
(bright enough to exhibit a crude spectrum), but only an 
upper limit of 0.004 IPC counts s ' with a 5 times longe r 
observation, nearly one year later jWalter & Kuhil |1984|) . 
Walter & Kuhi (1984) associated the X-ray emission detected 
from AATau as a quiescent level, and interpreted the non- 
detection in the framework o f the smothered coronae of 
TTSs (see IWalter & Kuhil 1198 Ih . where it was proposed that 
mass ejection could increase enough the absorbing column 
density to smother th e coronal (quiescent) X-ray emission. 
iNeuhauser et alJ (1 19951) reported 10 years later a RASS de- 
tection with 0.014 + 0.006 PSPC counts s->0 AATau was also 
detected during the PSPC pointed observation with 0.022 + 
0.005 PSPC counts s ' (see in the WGA catalogue of ROSAT 
poin t sources, the sourc e 1 WGA J0434.9-H2428, located 37' off'- 
axis; IWhite et al.l[T996l) . but (again) only during the shortest ex- 
posure. 

Further constraint can be derived on the X-ray variability of 
AA Tau by combining these multiple epoch observations with 
our better assessment of its quiescent level thanks to better 
spectra. For a plasma temperature ~23 MK and an unabsorbed 
X-ray luminosity in the energy band from 0.5 to 8.0 keV of 
~ 0.5 X 10""' ergs"' (equivalent to an unabsorbed X-ray flux 
of 2.1 X lO^'-'ergs"' cm"^), absorbed by a column density of 
~ 10^^ cm"^, we compute with PIMM^Ithat Einstein/IPC (0.2- 
4.5 keV) and i?OSAT/PSPC (0. 12-2.48 keV) would both observe 
only ~0.002 counts s"' . This low count rate is consistent with the 



^ iReale et ( 119971) use (, the slope of the flare decay in a log-log 
diagram of the plasma temperature vs. the squared-root of the emis- 
sion measure (a proxy of the plasma density), as diagnostic to assess 
the level of sustained heating in the analysis of stellar flares. For the 
XMM-Newton energy coverture, the corresponding range for f values 
are 0.4 and 1.9, for strongly sustained heating and freely decaying loop, 
respectively. 

Note t hat SR reported a RASS rate 10 times higher than the one 
reported bv lNeuhauser et all d 19951) . 
' Available at http : //asc ■ harvard . edu/toolki t/pirtims . j sp | . 
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IPC upper limit0 and we found that it is twice lower than the 
background level (inside a 1.5'-radius circle) at the location of 
AA Tau in the second PSPC pointed observation. Therefore, the 
previous non-detections in X-rays are consistent with the quies- 
cent level observed with XMM-Newton. 

We conclude that the previous X-ray detections with 
Einstein, the RASS, and ROSAT/PSPC, were made during high 
levels of activity, where AATau was about 18, 8, and 13 times, 
respectively, above its quiescent level. 

2.5. XMM-Newton optical/UV monitor light curves 
2.5.1. UV photometry 

The OM was operated in the imaging mode default, which uses 5 
imaging windows plus a small (1.7'xl.7') central imaging win- 
dow. The imaging mode default consists of a sequence of 5 expo- 
sures where one of the 5 imaging windows covers a large fraction 
of the OM field-of-view (17'xl7') with 1" x 1" spatial resolu- 
tion, while the small central imaging window ensures a contin- 
uous monitoring of the target at the center of the field-of-view 
with (y.'5 X (X'S spatial resolution (for an illustration of this ex- 
posure sequence and a light curve obtained with the small cen- 
tral window s ee Fig. 85 of the XMM-Newton Users' Handbook 
and iGrosso et al.l i2007 . respectively). We requested the mini- 
mum available exposure time in imaging mode default (800 s) 
to monitor any change in the UV photometry (UVW2 filter at 
206 nm) with a ~15 min sampling rate. 

We run the OM imaging mode pipeline. We made our own 
IDL program to plot the light curves of AA Tau from the source 
Hsts of the small centi-al window (*OM*SWSRLI0000.FIT), and 
the large central window (*OM*SWSRLI1000.FIT). The latter 
provides only 1 data point per OM exposure sequence (about 
5 X 800 s), which is obtained simultaneously with the first (out 
of 5) data point of the small central window. In a few exposures, 
AA Tau is missing in the observation source list of the small cen- 
tral window, but a visual inspection of the corresponding images 
confirms the detection. Therefore, we complete the light curve 
by doing aperture photometry. 

The top panel of Fig.[T]shows the UV light curves of AA Tau. 
We note that the UV light curves of AA Tau reported in SR were 
limited to the UV photometry obtained with the large OM central 
imaging window (see for comparison the grey hourglasses in the 
top panel of our Fig. [1] and their Fig. 1). The continuous moni- 
toring with the central window is crucial to determine accurately 
the UV variations on hour timescale. We find that AATau is 
variable in UV by an amount of a few 0. 1 mag on a few hours 
timescale, and up to 1 mag on a week timescale. There are no 
correlation between the UV and X-ray variations. 
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Fig. 5. Weekly and daily UV variability of AA Tau observed with 
the OM. Top panel: the dotted and dashed-dotted line show the 
weekly modulation of the UV flux attributed to the eclipse pe- 
riod (8.22 days), and the overall modulation of the UV flux. The 
horizontal line shows the average flux. Bottom panel: the dashed 
line shows the daily modulation of the UV folded in phase after 
subtraction of the weekly modulation. Symbols are the same in 
both panels. The formula of the overall fit (dashed-dotted Une) 
is given by Eq. \T\ 



reproduced with this simple model including only weekly vari- 
ations. A Lomb periodogram of the residuals suggests the pres- 
ence of an additional (high-frequency) modulation with a period 
of 0.87 + 0.03 day. 

We then compute a grid of fits in the amplitude and phase 
parameter space to look for optimized values that maximize the 
peak of the periodogram residuals. This is equivalent to a simul- 
taneous fitting of both modulations with the long-term period is 
fixed. We find a better fit (x^ = 103.8 for 77 d.o.f.). Given the 
new and old values of and number of degrees of freedom, a 
F-test indicates a probability of -10"" for the null hypothesis; 
therefore, we conclude that it is reasonnable to add this extra 
modulation to improve the fit. Moreover, a Lomb periodogram 
of the residuals shows no other modulations, and a perfect sub- 
straction of the two modulations. The fit formula is given by the 
following equation: 



2.5.2. Modeling of the UV variability 

The UV excess observed in CTTSs is attributed to the accretion 
shocks at the base of the acc retion funnel flows (se e review on 
magnetospheric accretion by iBouvier et al.ir2007bl) . Therefore, 
we assume that the UV flux must be somehow modulated by the 
warped inner disk, and we perform a least squares fitting of the 
UV light curve using a cosine functi on with free amplitud e and 
phase, and period fixed to 8.22 days CBouvier et al.l l2007a^. The 
resulting fit (x^ = 206.5 for 80 d.o.f.) is not acceptable, mainly 
due to fast variations on day timescale that cannot be properly 



However, the RASS rate is not consistent with the IPC upper hmit. 
Note that SR argued the opposite. 



UV = 16.30 + 0.02 mag (1) 
+ (0.20 + 0.02) cos{27r[r - (17.2 + 0.1)]/8.22) 
+ (0.24 + 0.02) cos{27r[f - (14.75 + 0.02)]/(0.87 + 0.04)), 

where t is the day in February 2003, and the errors are given at 
the 68% confidence level. 

The top panel of Fig. |5] shows the UV light curve and the 
fit given by Eq. [T] The bottom panel shows the UV light curve 
folded in phase with the high-frequency period after the sub- 
straction of the long-term modulation. This light curve looks 
rather convincing. However, a longer UV observation with the 
OM is necessary to have a definitive confirmation of this high- 
frequency period, which would suggest a non-steady accretion. 
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Table 3. Journal of ground-based optical observations. 
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Fig. 6. Optical light curves of the OM guide stars obtained with 
the field acquisition exposure. The dashed lines and grey stripes 
indicate the raw count rates of guide stars and one-sigma vari- 
ations, respectively. The thick line and error bars shows the op- 
tical light curve of AA Tau. Labels indicate counterparts in the 
Guide Star Catalog (version 2.3.2). 



2.5.3. Optical variability from the field acquisition exposure 

We propose here to obtain extra informations on the optical pho- 
tometry of the target thanks to the OM Field Acquisition expo- 
sures (FAQs), where several stars are also detected in optical in 
the OM field of view, and are then available as comparison stars 
for variability study. 

The FAQ is a short exposure (10 s) V-filter image always 
taken at the start of each science observation (i.e., at the begin- 
ning of the MOS observation) to allow proper identification of 
guide stars and to compensate small pointing errors. A threshold- 
ing is applied to the FAQ image by the OM on board software 
to identify guide stars. The derived offset is then applied to the 
OM science windows. The informations (positions in detector 
coordinates, count rates,...) on the guide stars are reported in the 
"OMS4O0OORFX" data files delivered with the Observation Data 
Files (ODFs). The use of the FAQ image, instead of scheduling 
an additional OM exposure in the imaging mode default with the 
V-filter, allowed us to make a twice longer observation with the 
UVW2 filter, by saving at least 5 x 800 s (without taking into 
account time overheads). 

In each observation, we recover sky coordinates from detec- 
tor coordinates, and identify the counterparts of guide stars in 
the Guide Star Catalog (version 2.3.2). Then, we build a light 
curve for the 13 guide stars (including AATau). Fig. |6] shows 
the variations of the raw count rates in the optical of the guide 
stars during our campaign. All guide stars, except AA Tau, ex- 
hibited small (lower than 2cr) relative variations in the optical. 
By contrast, AA Tau exhibited large relative variations, in par- 
ticular, during observations #1, #2, #5, and #6 where a large dip 
is visible. We conclude that these 4 XMM-Newton observations 
were likely made during the optical eclipses of AA Tau. 
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3. UV and X-ray variabilities versus rotational phase 

3.1. Ground-based optical photometry 

The Journal of ground-based optical observations is given in 
Table [3] Observations were carried out from three sites over a 
time span of two weeks, covering our XMM-Newton observa- 
tions in February 2003, using either CCD detectors or a pho- 
tomultiplier tube (Mt Maidanak). Measurements were obtained 
in the V filter. Differential photometry was performed on CCD 
images and absolute photometry from photomultiplier observa- 
tions, with an accuracy of the order of 0.01 mag. Somewhat 
larger systematic errors (<0.05 mag) might result from the rel- 
ative calibration of the photometry betwe en site s . All data re- 
duction procedures can be found in .Bouvier et al. I (l2003l) . 

The top panel of Fig. |7] shows the ground-based optical 
light curve; for comparison purpose, the other panels show 
the UV and X-ray light curves, and the plasma parameters. 
Unfortunately, due to bad weather conditions in Europe, none 
ground telescope was able to obtain simultaneous optical obser- 
vations with XMM-Newton. In the most favorable cases, ground 
optical photometry was obtained only 4.3 h, 3.1 h, 2.8 h, and 
2.0 h before the beginning of the MOS observations #l-#4, re- 
spectively; and only 3.3 h after the end the last MOS observa- 
tion. Despite the limited time sampling, a dip, likely associated 
with a primary eclipse of AATau, is however detected at the 
beginning of the monitoring campaign. This supports what we 
deduced previously in Sect. 12.5.31 from the optical FAQ images, 
that XMM-Newton observations #1 and #2 were obtained dur- 
ing the eclipse. Consequently, the FAQ images can also be used 
safely to identify an eclipse during XMM-Newton observations 
#5 and #6, located within the Feb. 23-27, 2003 gap of the ground 
observation. 



3.2. Folded light curves and eclipse phases 

To obtain a better determination of the dates of the optical 
eclipses, we compare our photometry with the on e obtained just 
6 months apart, from Aug. 27 to Oct. 24, 2003 (iBouvier et al.l 
i2007al) . The top panel of Fig. [8] shows both photometry data 
set f olded together in pha se using the 8.22-day rotational pe- 
riod dBouvier et al.ll20()7ah . i.e., we introduced no phase differ- 
ence between the two data set. The (arbitrary) common origin of 
phase was ch osen to match t he hgh t curve of season 2003 shown 
in Fig. 15 of iBouvier et al.l (l2007ah . taking the eclipse center at 
phase 0.5. 

The behaviour of our photometry is consistent with the shape 
of the light curve observed during the season 2003. We can ex- 
clude any phase offset larger than 0.05 between the two epochs. 
A deep primary eclipse (AV ~ 1.5 mag) is well detected at the 
beginning of our monitoring campaign. There is also some ev- 
idence for a shallow secondary eclipse detected only at the end 
of our monitoring campaign, after the last XMM-Newton ob- 
servation. The differences in brightness of about 0.5 and 1 mag 
between the Feb. and Aug. 2003 data, observed at the begin- 
ning and the end of the primary eclipse, respectively, can be 
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Fig. 7. Light curves and plasma parameters of AA Tau. The panels show from top to bottom: the optical, UV, and X-ray light curves 
(labels indicate observation numbers); the variations of the plasma emission measure(s) and temperature(s), and the photoelectric 
absorption (Table [2li . The top horizontal ax is indicates the corresponding phase for the rotation period of 8.22 days. The (arbitrary) 
phase origin is from lBouvier et all (l2007ah . 
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explained with a longer duration of the eclipse in Feb. 2003. 
The brightening on Feb. 16, 2003, close to the center of the pri- 
mary eclipse, was observed with a high sampling rate (~ 2 min), 
and exhibited only a fast decay. Therefore, this event is similar 
to the transient brightening events u sually observed in the faint 
state of AA Tau (iBouvier et al.ll999l) . The primary eclipses were 
centered on Feb. 15.5 and Feb. 23.5, 2003, and covered phases 
ranging from ~0.3 to ~0.7. Our ground-based monitoring con- 
firms that XMM-Newton observations #1, #2, #5, and #6 were 
made during primary eclipses. In particular, observations #2 and 
#6 are secured close to the center of the primary eclipse. We can 
exclude a secondary eclipse at the start of observations #4 and #8 
thanks to the FAQ, and moreover the corresponding X-ray light 
curves show no decay. Therefore, we conclude that the shallow 
secondary eclipse likely started after the end of our observations. 

The UV flux is the highest when the primary eclipse starts 
and the lowest towards the end of it. Indeed, the lowest UV flux 
was observed during the XMM-Newton observation #3 at the 
end of the egress phase, i.e., outside the primary eclipse. Our 
model of the UV flux variations helps to disentangle weekly 
modulation (produced by the warped disk) and daily variation. 
Eq.[T]indicates that the weekly modulation was at minimum on 
Feb. 17.2, 2003 (see also the top panel fo Fig. |5]), which cor- 
responds to a phase delay of about 0.2 compared to the opti- 
cal eclipse. Therefore, the warped disk produces a maximum 
of obscuration of the UV flux at the end of the optical eclipse. 
Consequently, the true maximum of the UV flux occurred around 
phase 0.2, well before the start of the eclipse. We didn't sur- 
vey this time interval with XMM-Newton, but we note that a 
brightening in the B-band of A A Tau was observed a round phase 
0.2 during the 1995 campaign (i Bouvier et al.lll999l) . The delay 
between the optical and UV eclipse, and the smaller depth of 
the UV eclipse (AL'yW2 ~ 0.40 mag) compared to the optical 
echpse (AV ~ 1 .5 mag), suggests a trailing accretion funnel flow, 
producing a strong absorption of the UV photons emitted at the 
accretion shock. 

No eclipses are detected in X-rays. The variable photoelec- 
tric absorption of the X-ray spectra is not correlated with the 
rotational phase; similar low and high values of the column den- 
sity are observed both during the eclipse and outside it. The 
maximum of the column density was observed at phase 0.6, 
close to the center of the optical eclipse, during the X-ray flare. 
However, increases of the column density were also reported 
during large stellar flares, and solar flares are sometimes associ- 
ated with coronal mass ejection (see the review on X-ray astron- 
omy of stellar coronae by Giidel 2004, and references therein). 
Therefore, we cannot rule out that the peak of column density 
is due to this energetic event. The gas column density on the 
line of sight p roduced by the warped disk is around 10^^ cm - 
dBouvier et al.. 1999), which is large enough to absorb all the X- 
rays emitted by AA Tau, even during a bright and hot flare. The 
lack of eclipses in X-rays indicates that X-ray emitting regions 
are located above the high-density disk warp at high latitudes. 

The online Fig.|9]shows the scattered plot of the average OM 
count rate versus the column density, where the count rate error 
include the observed variations of the UV flux (see Fig.[T]). The 
correlation coefficient of this sample is -0.72, suggesting a true 
correlation between the two physical parameters as argued by 
SR, who used smaller error bars for the count rates. However, we 
note that the data obtained outside the eclipse are located within 
the cloud of points, with error bars covering all the range of ob- 
served values of column density and count rate. Therefore, the 
column density variation cannot be attributed to the disk warp. 



Assuming a frequency of one X-ray flare per 6 50 ks, as ob- 
served in young solar-mass stars dWolk et alj|2005h . our 38.8 h- 
exposure should detect only 0.2 flare from AATau. The center 
of the optical eclipse is limited to about 0.2 in phase. Therefore, 
the combined probability to observe by chance an X-ray flare 
during the center of the transit of an accretion funnel flow is 
only 0.04. Moreover, the median flare level in young solar-mass 
stars is only 3.5 times the characteristic level ( Wolk et al. 2005), 
whereas we observed a flare with a larger amplitude. This sug- 
gests that this event is associated with the magnetic area corre- 
sponding to the base of the dipolar magnetic field line, which 
controls the accretion funnel flow. 



4. Discussion and conclusions 

The high throughput of XMM-Newton allows us to obtain spec- 
tra of AA Tau at each phase of its activity. By contrast to previ- 
ous i?OSAT observations, where the measurement for this source 
was limited to count rate and hardness ratios in the soft X-ray en- 
ergy bancQ, the plasma parameters can be derived from spectral 
fitting, which provides, in particular, a better measurement of the 
column density and the X-ray luminosity. 

We find that the column density, derived from the photo- 
electric absorption of the X-rays emitted by the active corona 
of AA Tau, varies from A^h ~ 1 -0 x 10^^ cm"^ to 1 .8 x 10^^ cm"^. 
However, the optical extin ction of AA Tau by the dust is very 
low with Ay = 0.78 mag (IBouvier et alj|1999l) . which can be 
converted to A^h ~ (1.2 + 0.1) x 10^' cm~ ^, using the relatio n 
Nii/A, = (5.6 + 0.4) X lO^^cm-^mag-^ (Vuong et al. '2003b, 
combined with the extinction law of Rieke & Lebofskv ( 1985]), 
A] - 0.282 X Ay (for their adopted Ry value of 3.1). We con- 
clude that there is an excess of column density, varying from 
0.9 X 10^^ to 1.7 X 10^^ cm"^. These values are consistent with 
the one reported by SR. To explain this excess of column density, 
SR introduced an additional absorption in a disk wind (with no 
dust, and cool enough to avoid producing any soft X-ray emis- 
sion), or a peculiar dust grain distribution with Ry ~ 0.4 to rec- 
oncile the observed extinction and X-ray absorption. However, 
we note that three-dimensional MHD simulations of ^isk accre- 
tion to a rotating magnetized star (e.g., Ro mano va et alj|2003l) 
show that the stellar magnetosphere is far to be an empty cav- 
ity. Matter accretes mainly through narrow funnel loaded with 
high-density material (~10'^cm"^), which are surrounded by 
lower density funnel flows that blanket nearly the whole mag- 
netosphere. Therefore, we propose to identify this excess of gas 
with this lower density funnel flows filling the stellar magne- 
tosphere. Dividing the excess column density by the width of 
the magnetosphere of AATau (about 7.8 stellar radii) leads to a 
density of a few 10'" cm"-', which is compatible with this inter- 
pretation. Moreover, the multiple spirals visible in the simulated 
accretion flows should help to produce a variable column den- 
sity. This low density gas, located below the radius of dust sub- 
limation (close to the inner accretion disk for AA Tau), is then 
dust free. 



^ To convert the RO SAT count rate of AA Tau to X-ray luminosity, 
iNeuhauser et al.l ( 11995 ') estimate an energy conversion factor from the 
visual extinction (converted to foreground hydrogen column density) 
and the second hardness ratio (sensitive to the plasma temperature). 
This leads to an X-ray luminosity of 0.4 x 10^" ergs"' (see Table 1 of 
Johns-Rrull 2007), which is likely underestimated by a factor of ten (see 
our simulations of count rates in Sect. l2.4t . because the column density 
of AATau (see Table[2l and di scussion below) i s about 10 times larger 
than the one assumed bv iNeuhauser et alJ (119951) . 
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Accreting TTSs show in average a Lx/^boi ratio 2.5 times 
smaller than in non-accreting TTSs (see for COUP and XEST re- 
sults, Preibisch et al. 2005 and Briggs et al. 2007, respectively). 
This is generally interpreted as a direct or indirect consequence 
of the accretio n proces s (see discussion inlPreibisch et al.ll2005t 
iTelleschi et al . 2007; Gregorv et al.ll2007h . The median value of 
log(Lx/i'boi) for AA Tau is -3.5; for comparison this value is in- 
termediate between -3.7 and -3.3, the me dian values found for 
accretors and non- accretors, respectively (IPreibisch et al.ll2005t 
iBriggs et al.l l2007l) . The peculiar orientation of AATau may 
help to overcome pa rtly the extinction by the accretion funnels 
jGregorv et al.ll2007l). 

Recently, " iJohns-KrullI (l2007h reported new magnetic field 
measurements for CTTSs, based on Zeeman broading of pho- 
tospheric absorption lines in the near-IR, showing that the ob- 
served mean magnetic field is in all cases grea ter than the 
field predicted by pressure equipartition arguments. Johns-Krull 
(12007) suggests that the very strong fields decrease on these stars 
the efficiency with which convective gas motions in the photo- 
sphere c an tangle magnetic loops in the corona. For AATau, 
iJohns-Kr ull (2007) reports a value of 2.78 kG for the mean mag- 
netic field, that is 2.7 times larger than the field strength at 
pressure equipartition; and predicts from solar scaling an (un- 
observed) X-ray luminosity of about 4 x 10^" ergs"' (for com- 
parison this value is 10 times greater than the minimum level 
that we observed). However, our observation shows that a bright 
X-ray flare can occur during the transit of an accretion funnel 
flow, and locate the active X-ray area likely close to the foot of 
the accretion funnel flow. We speculate that a magnetic interac- 
tion exists between the free falling accretion flow and the strong 
magnetic field of the stellar corona, which may trigger magnetic 
reconnections and give rise to bright X-ray flares, which boost 
the X-ray emission. 

This campaign of observations of AATau with XMM- 
Newton shows that X-ray spectroscopy with CCD provides a 
unique tool to probe the circumstellar dust-free gas in this ob- 
ject, and that some magnetic flares are likely associated with the 
accretion process. Longer coordinated optical and X-ray/UV ob- 
servations of AA Tau are still needed to obtain continuous light 
curves, which is crucial to confirm the flaring behaviour of ac- 
tive regions associated with the accretion funnel flow, and daily 
modulation of the UV flux. Simultaneous Zeeman-Doppler im- 
ages would help to derive a surface magnetogram; this map of 
the magnetic active regions on the stellar photosphere, combined 
with the X-ray light curves, would then allow to build a self- 
consistent three-dimensional model of the corona of AA Tau. 

Taking into account that the duration of the visibility window 
for AATau is about 130ksec per XMM-Newton orbit (2 days), 
a full monitoring of two optical eclipse periods of AA Tau (8.22 
days) with XMM-Newton is equivalent to an effective exposure 
of 1.1 million seconds. This project has the typical duration of 
the large programs which are anticipated with XMM-Newton for 
the next decade. 
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Appendix A: Elemental abundances of the coronal 
plasma and the absorbing material 

The typical values of plasma element abundances observed 
in the coronae of young stars with fine X-ray spectroscopy 
dGudel et al.ll2007ah . that we use for the vapec coronal plasma 
model in our X-ray spectral fitting with XSPEC, are given in 
Col. (2)-(4) of Table ED 



Table A.l. Elemental abundances of the coronal plasma and the 
absorbing material. 



The column density of the absorbing material located on the 
line of sight is estimated from spectral fitting using a photoelec- 
tric absorption model, that uses photoionization cross sections, 
and solar abundances for the material composition. The wabs 
photoelectric absorption model use photoionization cross sec- 
tions of Morrison & McCammon ( 1983), and (old) solar abun- 
dances of lAnders & Grevessel(ll989() . Significant revisions of the 
solar abundances have been made recently, as a result of the ap- 
plication of a time-dependent, 3D hydrodynamical model of the 
solar atmosphere, instead of ID hydrostatic models. This has de- 
creased the metal abundances, in particular of carbon and oxy- 
gen, which are the main contributor to the photoionization cross 
section above 0.3 and 0.6 keV, respectively. Consequently, the 
absolute value of the column density is dependent of the adopted 
photoelectric absorption model. Using updated solar abundances 
is then crucial when the absolute value of the column density 
is needed (e.g.. IVuong etaPIIOOl) . Col. (5)-(7) of Table \AJ\ 
give the recent comp ilation of elemental solar abundances by 
lAsplund et al.l (2005), where the de crease of metal by compari- 
son with lAnders & Grevess^ (Il989l) is indicated in the last col- 
umn. 

We test the impact of these udpated s olar abundances o n 
our fitting by replacing wabs by tbvarabs (IWilms et al.ll2000h . 
which allows to input new abundances for the absorbing ma- 
terial. Moreover, tbvarabs uses also updated photoionization 
cross sections. We find nearly identical values of temperature 
and emission measure, however, as anticipated the column den- 
sity value is increased. Fig. lA. II shows that column density val- 
ues obtained with wabs are underestimated by about 50%. 



Coronal plasma" 



Absorbing material'' 



El A(E1) 
(1) (2) 



«(El)/n(H) 
(3) 



angr 
(4) 



A(E1) 
(5) 



n(El)/«(H) 
(6) 



angr 
(7) 



He 
C 
N 
O 
Ne 
Na 
Mg 
Al 
Si 
S 
CI 
Ar 
Ca 
Cr 
Fe 
Co 
Ni 



10.99 
8.21 
7.95 
8.56 
8.01 



7.00 
6.17 
7.04 
6.83 

6.30 
5.65 

6.96 

5.54 



9.77E-02 
1.63E-04 
8.83E-05 
3.63E-04 
1.02E-04 

9.99E-06 
1.47E-06 
l.lOE-05 
6.76E-06 

2.00E-06 
4.47E-07 

9.13E-06 

3.47E-07 



1.000 
0.450 
0.788 
0.426 
0.832 



0.263 
0.500 
0.309 
0.417 



0.550 
0.195 



0.195 
0.195 



10.93 
8.39 
7.78 
8.66 
7.84 
6.17 
7.53 
6.37 
7.51 
7.14 
5.50 
6.18 
6.31 
5.64 
7.45 
4.92 
6.23 



8.51E-02 
2.45E-04 
6.03E-05 
4.57E-04 
6.92E-05 
1.48E-06 
3.39E-05 
2.34E-06 
3.24E-05 
1.38E-05 
3.16E-07 
1.51E-06 
2.04E-06 
4.37E-07 
2.82E-05 
8.32E-08 
1.70E-06 



0.871 
0.676 
0.538 
0.537 
0.562 
0.691 
0.892 
0.795 
0.912 
0.852 
1.682 
0.417 
0.892 
0.902 
0.602 
0.967 
0.954 



Notes: Col. (2) and (5) give the element abundances on the logarithmic 
astronomical scale, where the numbers of hydrogen atoms are set to 
A(H) = log)i(H) = 12. The numbers of element atoms normalized to 
the number of hydrogen ato ms are given in Col . (3) a nd (6); Col. (4) 
and (7) compare this ratio to I Anders & Grevessd ( Il989l) 's photospheric 
abundances. 

" Elemental abundances o bserved in the coron ae of young stars with 
fine X-ray spectroscopy iGildel et al.ll2007a|) used i n vap ec. 
Elemental solar abundances of " lAsplund et all ( |2005|) used in 
tbvarabs. 
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Fig. A.l. Comparison of the column dens ity values obtained 
from spectral fitting w hen using the old ( Anders & Grevessd 
Il989l) and the updated jAsplund et al.l 12005 ) solar abundances. 
The dashed line shows the mean average between the two values 
of column density. 
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Fig. 2. A subset of X-ray flares from the Chandra Orion Ultradeep Project with peak amplitude and duration larger than the one 
observed in AA Tau. The left panels show COUP light curves (see Getman et al. 2005b), where blue and red segments indicate 
the minimum (BBmin) and maximum (BBmax) levels obtained from Bayesian block analysis dScargldl 19981) . respectively. The blue 
dotted line show the minimum level. The peak amplitudes are also given. The dashed green lines start 1 .7 days after the end of the 
maximum Bayesian block, and indicate the level (2.5 times above the minimum level), that was observed in the third observation of 
AA Tau. The right panels are an enlargement of the light curve around the Bayesian block showing the peak level and duration. Dots 
mark the arrival times of individual X-ray photons with their corresponding energies given on the right-hand axis. Large vertical 
gray stripes indicate the five passages of Chandra through the van Allen belts where ACIS was taken out of the focal plane and thus 
was not observing Orion. 
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Fig. 3. EPIC pn (black), MOSl (red), MOS2 (green) spectra of AATau plotted with the same scale. The lines show our best fits 
using one-temperature plasma combined with photoelectric absorption for observations #1 to #8 (Table|2]l from left to right and top 
to bottom. The residuals are plotted in terms of sigmas with error bars of size one. 
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Fig. 9. Average OM count rate versus column density. Labels indicate XMM-Newton observations. The crosses mark observa- 
tions obtained outside the optical eclipse. The dotted hourglass shows the column density of observation #2 when using only 
one-temperature plasma. 



